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Abstract. We present Keck R and I band images of the discovered: 1RXS J170849. 0—400910 (Sugizaki et al 
field of the anomalous X-ray pulsar IE 2259+58.6. We 1997|); IE 1841-045 in the supernova remnant Kes 73 QVa- 



derive an improved X-ray position from archival ROSAT sisht fc Gotthelf 1997); and AX J1845. 0-0300 (Torii et al 



HRI observations by correcting for systematic (boresight) 1998; Gotthelf fc Vasisht 1998) in the supernova remnant 



errors. Within the corresponding error circle, no object 
is found on the Keck images, down to limiting magni- 
tudes R = 25.7 and I = 24.3. We discuss the constraints 
imposed by these limits, and conclude that it is unlikely 
that IE 2259+58.6 is powered by accretion from a disk, 



G29. 6+0.1 JGacnslcr ct al. lgff ). Indeed, perhaps even 
the progeny of these systems ha s been identified, in RX 
J0720.4-3125 ( |Habcrl ct al. 1997j ). The latter object is the 
only one for which the optical counterpart has been iden- 
tified flMotch fc Habcrl 1998]; |Kulkarni fc van Kerkwijk 



irrespe :tive of whether it is in a binary or not, unless the 
binary is extremely compact. 

Key words: pulsars: individual: IE 2259+58.6 - accre- 
tion, accretion disks - Stars: neutron - X-rays: stars 



1998|) 



A variety of models for AXPs has been put forward, us- 
ing both accretion and intrinsic mechanisms to explain the 
X-ray luminosity. For the accretion scenarios, the source of 
matter has been suggested to be a very low mass compan- 



ion (Mcrcghctti & Stella 1995) or a circumstellar disk com 



1. Introduction 



posed of debris, either from a common envelope phase of 
a high mass X-ray binary ( van Paradijs et al. 1995]; phosh 



There is a growing number of X-ra}r pulsars for which it 

is not clear why they emit X-rays. Their rotational energy 
loss appears insufficient to power the X-ray luminosity and 
there is no sign of a companion from which they could 
accrete matter. These pulsars are called the "anomalous 
X-ray pulsars" (AXP). 



ct al. 1997) or from a supernova explosion (Corbet ct al 
1995; Chatterjee et al. 1999). In scenarios in which the 



X-ray emission is intrinsic, the sources are either hot, mas- 
sive, rapidly spinning white dwarfs, formed in the merger 



of a double white-dwarf binary (Paczyhski 199C), or so- 
called "magnetars," isolated neutro n stars with very high 
magnetic field strengths (<; 10 14 G; Thompson fc Duncan 



Th e AXP diff er from other X ray puloaro in the fol 
lowing properties (|Hcllicr 1994 tVIcrcghctti fc Stella 1995| . 
van Paradijs et al. 1995 ): (i) their spin periods are within 
a small range, between 5 and 9 seconds, and their spin- 
down rates are nearly constant; (ii) timing analysis gives 
no evidence of orbital motion; (iii) their X-ray spectra are 
rather soft; (iv) their X-ray light curves show less vari- 
ability; (v) their inferred X-ray luminosities are rather low 
(~ 10 35 ergs -1 ); (vi) they are located close to the Galactic 
plane; (vii) some appear associated with supernova rem- 
nants. 

Members of the group of AXP include 4U 0142+61, 
IE 2259+58.6, IE 1048.1- 5937 and RX J1838.4-0301 
(Mereghetti & Stella 1995 and references therein). Re- 



1996). In the latter case, AXP might be related to soft 



cently, three X-ray pulsars with similar properties were 



gamma-ray repeaters]]. 

In this paper, we present a deep search for the opti- 
cal counterpart to the AXP IE 2259+58.6. In Sect. |, we 
briefly review what is known about this source. We de- 
rive a refined position for IE 2259+58.6 in Sect. 0, and 
in Sect. | we present our optical observations. We discuss 
the results in Sect. |[ 

2. IE 2259+58.6 

IE 2259+58.6 is located inside the supernova remnant 
CTB 109 (G109.1-1.0). This remnant is semi-circular, 
with neither X-ray nor radio emission from its western 
side (see, e.g., Rho fc Pctrc 1997), which may be related 
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1 The proposed X-ray co unterparts to SGRs also have spin 
periods in the 5-9 s range (Kouveliotou et al. 1998, 1999). 
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to the presence of a molecular cloud. A lobe of enhanced 
X-ray emission ('jetlike' feature) is present in between the 
eastern shell and the pulsar, but likely this does not reflect 
physical interaction. 

IE 2259+58.6 is an X-ray pulsar with a period of 6.98 s. 
It is observed to be spinning down with a period derivative 



of 5.5 x lO-^ss" 1 . Baykal & Swank ( 199J ) found small 



but significant deviations from constant spin-down and 
concluded th at acc retion must be the source of energy (but 
see Melatos 1999 for a different explanation). An upper 
limit ax sin i < 0.03 lt-s has recently been obtained with 
the RXTE ( [Mereghetti et al. 199g| ). Assuming a Roche- 
lobe filling companion and an orbital inclination i > 30°, 
this restricts the companion mass to ^O.IMq for a main- 
sequence star, and to ^0.8Mq for a helium-burning star. 
No constraint can be set for a white-dwarf companion. 

IE 2259+58.6 has a very soft X-ray spectrum, best 
described by a combination of a power-law with photon 
index 4.0+0.1 and a blackbody with kT = 0.43+0.02 keV, 
absorbed by a column density TVh = 8.5 x 10 21 cm -2 (Rho 



& Peti je 1997| using ASCA, BBXR T and ROSAT PSPC 
data; see also Corbet et al. 1995). A possible problem 



is contamination of the spectrum by extended emission 



around the pulsar (Rho & Petre 1997). Parmar et al. 
( |1998| ) tried to account for this by subtracting the ob- 
served spectrum of the supernova remnant multiplied with 
a scale factor. They find that the residual spectrum can 
be fitted with a simple power law with an absorption col- 
umn density of Nh = (2.18 ± 0.07) x 10 22 cm" 2 . It should 
be noted, however, that the BeppoSAX LECS and MECS 
instruments used by these authors have rath er po or spa- 
tial resolution. Furthermore, Rho & Petre ( 1997 ) found 
that the spectral shape of the extended emission near IE 
2259+58.6 differs from that of the remnant. 

No radio counterpart to IE 2259+58.6 has been found 



down to an upper limit of 50 /iJy (Coe et al. 1994). So far, 
no counterpart at optical and infrared wavelengths could 
be identified ( pavies fc Coe 1991| |Coe fc Pightling 1998| ). 



These authors obtained limiting magnitudes of B — 25, V 
= 24, R = 24.5, I = 23, H = 18.5, J = 19.6 and K = 18.4, 
although they could not exclude several candidates due 
to the fact that the X-ray position of IE 2259+58.6 was 
not known accurately enough. Note that their astrom etry 
differs from ours and that of Fahlman et al. (11982D see 
Sect. | and Davies & Coe |l99l| ). 



2.1. Dv,LuiiLti and itddnniny Lu IE 2259+58.0 



All distance estimates to IE 2259+58.6 use distances to 
either CTB 109 or the molecular cloud, i.e., it is assumed 
that the pulsar is at the same distance as CTB 109 and 
that the molecular cloud is either at the same distance as 
well or in front (and causes the absence of emission from 
the western part of the CTB 109). The distance to CTB 
109 was estimated using the s urface brightness di amete r 
relation by Sofue et al. ( 1983] ) and Hughes et al. ( 1984 ), 



at 4.1 and 5.6 kpc, respectively. A lower limit to the dis- 
tance of d ^ 6+ 1 kpc was obtained by Strom (priv.comm) 
from hydrogen 21-cm line absorption measurements to- 



wards the remnant (Braun & Strom 1986) 



For the molecular cloud, a kinematic distance of 5 ± 
1 kpc can be derived from the molecular line measure- 



ments of Kahane et al. (1985) (here, we used the linear 
fit to the galactic rotation curve of Fich et al. 1989, with 
rotation constants 9o = 220 kms -1 and i?o = 8.5kpc). 
Distance estimates have also been made for several Hn 
regions associated with the molecular cloud. Spectropho- 
tomctric parallaxes for the exciting stars were obtained by 



Crampton et al. (1975), leading to distances of 5.5, 5.6, 3.6 



and 4.0 kpc for the Hn regions S148, S149, S152 and S153, 
respectively. Furthermore, a lower limit of d <; 6 ± 1 kpc 
towards S152 was obtained by Strom (priv.comm) using 
hydrogen 21-cm line absorption. From the above, assum- 
ing that the fractional errors in the spectrophotometric 
parallaxes are about 25% (corresponding to uncertainties 
of about 0.5 mag in the absolute magnitudes inferred from 
the spectral types), we conclude that a distance to the 
molecular cloud of 5 ± 1 kpc is most likely, and, therefore, 
that IE 2259+58.6 is at d ^ 5 ± 1 kpc. 

The reddening to IE 2259+58.6 can be estimated from 
the X-ray hydrogen column density using the relation 
N K = (1.79 + 0.03) x 10 21 Ay cm" 2 , where Ay is the visual 



extinction in magnitudes ( Prcdchl fc Schmitt 1995 ) . From 
the fits to the X-ray spectrum of Rho & Petre ( 1997 ), 
values of A*h between 0.8 and 1.2 x 10 22 cm~ 2 are found, 
with a best fit value of 0.85 x 10 22 cm~ 2 . Thus, one infers 
Ay = 4.5 to 6.7 mag, with a best-fit value of 4.7 mag. For 
N K = 2.18 x 10 22 cm~ 2 , as found by Parmar et al. ( [1998D , 
one would infer Ay — 12.2 mag. As noted above, however, 
we feel this measurement is rather uncertain. 

The reddening can be estimated indirectly using CTB 
109. Fesen & Hurford (1995) measured Ha and H/3 rela- 
tive line intensities in optical spectra of 5 filaments, and 
found Ay = 2.5 to 3.7 mag, with substantial variations 
across the remnant. These values appear somewhat lower 
than those derived above from the X-ray absorption for 
IE 2259+58.6, and thus one may wonder about the asso- 
ciation, or about possible intrinsic absorption around IE 
2259+58.6. Since the reddening to the remnant is seen to 
vary, however, a somewhat higher value for IE 2259+58.6 
does not seem unlikely. Furthermore, from the X-ray emis- 
sion of CTB 109, column densities corresponding more 



close ly to that of IE 2259+58.6 are inferred (Rho fc Petre 



1997) 



In the following, we will adopt a reddening of Ay = 
4.7 mag and a distance of 6 kpc, unless noted otherwise. 

3. Improved X-ray position 

IE 2259+58.6 has been in the field of view of two ob- 
servations with the ROSAT high resolution imager (HRI; 
David et al. 1999). The first was in 1992 from January 8, 
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07:16 UT, to January 10, 18:59 UT, and was centered on 
IE 2259+58.6. The second was pointed at the lobe to the 
East of the pulsar and was from June 25, 07:46 UT, to 
June 27, 08:26 UT. In both observations, the semicircular 
shell, the lobe, the pulsar and the extended emission in 
its immediate surroundings, as well as three other point 
sources can be seen. We identify the latter with sources 1, 



Table 1 . X-ray and optical positions used for the astrom- 
ctry. 



2 and 3 in Table 4 of Rho & Petre ( [19971 ), and will refer 
to these as RP1, RP2 and RP3 hereafter. 

The data were analysed using the EXSAS software 
package (1998 April version) in the following way. First, 
we corrected for the small error in HRI plate scale found 
by Hasinger et al. ( 1998j ), by multiplying all detector posi- 
tions (relative to the detector centre) with a factor 0.9972. 
Next, we ran the standard source detection program, 
which creates a preliminary source list using a sliding cell 
algorithm. This source list is fed into a maximum likely- 
hood detection algorithm (Cruddace et al. 1988), which 
yields precise positions (see Zimmermann et al. 1998). In 
order to prevent trouble with the extended emission from 
the remnant, we used two masks, one to the East of CTB 
109 containing RP3 and one to the West containing RP1, 
RP2 and IE 2259+58.6. Our final positions are listed in 
table |l|; the 1-a errors are internal, i.e., they do not yet 
include pointing uncertainties. We verified that different 
choices of mask images gave consistent results. 

Comparing the positions derived for the two observa- 
tions (Table |l|), it is clear that they are not consistent. 
This reflects the uncertainties in the attitude solutions 
used to transform detector positions to sky positions (see, 
e.g., David et al. 1999). One can correct for this if one has 
independent positions for some of the sources, for instance 
from optical counterparts. Below, we will try to do this. 

Many X-ray sources can b e iden tified with bright stars; 
for instance, Motch et al. ( |l998|) find that 85% of the 
sources in a galactic plane region are bright stars. There- 



fore, we searched the USNO-A2.0 catalogue ( Monet et al 
1998()j3r bright stars (r usnoa2 < 13.0) near the positions of 
RP1, RP2 and RP3. We tentatively identify RP1 and RP3 
with stars USNO-A2 1425-14436845 and USNO-A2 1425- 
14515707, respectively. No bright star was found within 
30" of the position of RP2. We derived proper motions 
for the two stars by measuring positions on first and sec- 
ond generation images from the digitized sky survey, us- 
ing astrometry relative to the USNO-A2.0 catalogue^ (Ta- 
ble 0). For USNO-A2 1425-14515707, we found a signifi- 
cant proper motion of 0'.'022±0'.'004 yr _1 at position angle 
90?7± 10?2. 

The refined position for IE 2259+58.6 was derived as 
follows. First, the HRI positions were shifted to a common 



2 Note that both DSS1 and USNO-A2.0 positions were de- 
rived from the same plate and the results should therefore 
agree. The slight discrepancy in declination for RP3 is probably 
due to a second (fainter) star close to the proposed counterpart. 
This does however not affect our results. 
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RP1 - HRI 




















January 1992 


23 


00 


32.85 
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52 
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Tune 1992 


23 


00 


33.98 


58 


52 


49 
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2 
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1425-14436845 b 


23 


00 


33.396 


58 


52 


47 


16 





2 


DSS1 1953.830 


23 


00 


33.393 


58 


52 


47 


22 





2 


DSS2 1991.671 


23 


00 


33.377 


58 


52 


47 


07 
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RP3 - HRI 




















January 1992 


23 




i n a i 

iy.4i 


r, o 


4o 


v -> 

oo 


n 
9 


i 

1 


Q 
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June 1992 


23 


03 


20.30 


58 


45 


26 


4 


1 
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RP3 - counterpart 




















1425-14515707 b 


23 


03 


19.472 


58 


45 


28 


95 





2 


DSS1 1953.830 


23 


03 


19.471 


58 


45 


28 


43 





2 


DSS2 1991.671 


23 


03 


19.578 


58 


45 


28 


68 
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IE 2259+58.6 - HRI 




















January 1992 


23 


01 


08.09 


58 


52 


48 


3 





1 


June 1992 


23 


01 


08.99 


58 


52 


45 


3 





1 


IE 2259+58.6- corrected 


















Using RP1 and RP3 


23 


01 


08.44 


58 


52 


44 


1 





9 


Using RP1 only 


23 


01 


08.62 


58 


52 


43 


5 


1 


3 



a Uncertainty in each coordinate (in arcsec). 
b USNO-A2 identifier; mean epoch 1953.830. 



reference using the positions of IE 2259+58.6. Next, we 
determined weighted average offsets in right ascension and 
declination between the HRI and optical positions (includ- 
ing proper motion). For this purpose, we weigh the four 
individual offsets in right ascension and declination (RP1 
and RP3; two HRI observations) using the uncertainties 
given by the source detection algorithm (the errors in the 
optical positions are much smaller, ~0'/2). 

The differences with the weighted average correspond 
to x 2 = 11.5 for 6 degrees of freedom. This is not a good fit 
and suggests the presence of additional systematic errors. 
Following Hasinger et al. ( 1998 ), we added an error E sys 
in quadrature to the individual measurement errors such 
that Xr C d — 1- For this purpose, we require E sys — l'.'O. 
Including this additional uncertainty, we infer that our 
corrected position of IE 2259+58.6, listed in Table |l|, has 
an uncertainty of Of! 9 in each coordinate. This corresponds 
to a 95% confidence error radius of 2'.' 2. 



The value of -E S vs we 



by Hasinger et al. 



p98j ) 



find is twice the value used 
While these authors also in- 



cluded an 0'.'5 uncertainty in optical positions, we won- 
dered whether the somewhat large required value might 
be related to the rather large offset from the center of 
RP3. To get a feeling for the uncertainty, we also deter- 
mined a corrected position using RP1 only. The position 
is shifted slightly and has a slightly larger uncertainty of 
l'/3 in each coordinate (see Table Q). The 95% confidence 
error circle for this case has a radius of 3'.'3. 
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Table 2. Apparent magnitudes' 1 of selected stars in the 
field around IE 2259+58.6. Star designations are after 
Davies & Coe Jl99l] ). 



Star 




R 








I 










R 


- 1 


A 


19 


697 


± 





003 


18 


577 + 





005 


1 


120 


±0 


006 


B 


19 


265 


± 





003 


18 


217 + 





005 


1 


048 


±0 


006 


B' 


21 


206 


± 





010 


19 


987 + 





011 


1 


219 


±0 


015 


D 


20 


409 


± 





004 


19 


274 + 





005 


1 


135 


±0 


006 


E 


22 


537 


± 





016 


22 


545 + 





043 


-0 


008 


±0 


046 


F 


18 


902 


± 





004 


17 


904 + 





007 





998 


±0 


008 


H 


22 


398 


± 





014 


21 


219 + 





014 


1 


179 


±0 


020 


I 


23 


024 


± 





023 


21 


672 ± 





020 


1 


352 


±0 


030 


J 


2.3 


131 


± 





027 


21 


473 + 





019 


1 


658 


±0 


033 


K 


23 


676 


± 





037 


21 


297 + 





015 


2 


379 


±0 


040 


L 


23 


389 


± 





028 


21 


799 + 





028 


1 


590 


±0 


040 


M 


24 


557 


± 





086 


22 


832 ± 





068 


1 


745 


±0 


110 


N 


24 


848 


± 





110 


22 


119 ± 





031 


2 


729 


±0 


114 



a In addition to the listed uncertainties, there is a 0.1 mag 
uncertainty in the zero point; see text. 



We show both error circles superposed on Keck images 
in Fig. |]. Also indicated are the positions determined pre- 
viousl y by F ahlman et al. ( 1982 ; short dashed) and Hanson 
et al. ( |1988| ; long dashed) . Note that all these positions dif- 
fer significantly from that of Coe & Pighting (1998). We 
do not understand the reason for this discrepancy. 



4. Keck Observations 

Optical (R and I band) images of the field around IE 
2259+58.6 were obtained on 1994 October 30/31 and 1997 
January 7 using the Low-Resolution Imaging Spectrome- 
ter (LRIS, Okc et al. 1995 ) at the Cassegrain focus of the 
Keck I (1994) and Keck II (1997) telescopes. The 1994 
observations were plagued by clouds and high cirrus; the 
seeing was l'.'l. All 1994 images have a large gradient in 
the sky background. As a result, all but one of these im- 
ages are useless. In 1997, the weather was photometric, 
but the seeing mediocre at l'/2. All images were bias sub- 
tracted and flat fielded using dome flats. Unfortunately, 
it turned out that in the 1997 R-band images the posi- 
tion of IE 2259+58.6 was very close to a region of lower 
sensitivity in the CCD, which did not flat-field out. The 
resulting "smear" is just over stars K and N in two of the 
three images; it can also be seen in Fig. [l]). 

Astrometry was done using a short (300 s) R band im- 
age from 1997. We identified objects with stars in the 
USNO-A2.0 catalogue and obtain centroids for objects 
that were not saturated and appeared stellar (i.e., had 
a sufficiently small diameter). Next, we applied a bi-cubic 
distortion correction (J. Cohen, 1997, private comm.), and 
solved for zero-point position, scale, and position angle 
by comparison with the USNO-A2.0 positions. After re- 
jecting six outliers (probably misidentifications or large 
proper motion objects), the remaining 134 objects gave a 



good solution. The root-mean-square displacement from 
the catalogued positions is 0.26" in right ascension and 
0.25" in declination. 

For the pho tometry, we u sed the DAOPHOT II soft- 
ware package (stetson 1987). Instrumental magnitudes 



were obtained from two R band images (one from 1994 
and one from 1997, 1000s exposure each) and the average 
of two I band images (300s each, 1997 only). We cali- 
brated the photometry using the magnitudes determined 



by Davies & Coe (1991). The apparent magnitudes of se- 
lected objects as determined from the 1997 observations 
are listed in Table ||. The values derived for the 1994 R- 
band image are consistent with these, but have lower ac- 
curacy. Note that the uncertainties listed in the table are 
the formal la fitting errors. In addition to these, there 
are systematic unce rtaint ies related to the zero point. 
Since Davies & Coe ( |l99l| ) found a discrepancy by about 
0.15 mag in magnitudes obtained from two different expo- 
sures, the zero-point uncertainty is probably of the order 
of 0.1 mag. 

No object is found inside either of the HRI error circles 
derived above. We derived lower limits with the help of 
simulations. For these, we took the mean of the three 1997 
R band images and selected the inner region of the image 
that is void of stars. Next, we placed an artificial star 
with magnitude M at a random spot in this region and 
determined its magnitude in the same way as for the real 
objects (this also allowed us to verify the errors quoted in 
the table) . This was done ten times for a number of values 
of M, within a range 25.0-26.9 for R and 23.7-25.7 mag 
for I. We define the 3cr detection threshold to be that 
magnitude M for which the standard deviation equals 0.3 
mag. Taking account of the fact that we need about 3 
trials to cover the error circle, the derived 2a limits are 
R > 25.7 and / > 24.3 (for the error circle derived using 
RP1 only we need 7.5 trials and obtain R > 25.6 and 
/ > 24.2). 

5. Discussion 

We now examine the limits on the optical emission im- 
posed by the various models proposed for IE 2259+58.6. 
Contributions to the optical flux can come from the object 
itself, an accretion disk and a companion star. We will dis- 
cuss these, in the context of the proposed models, in the 
following subsections. 

5.1. An isolated object 



Paczyhski (1990) proposed that IE 2259+58.6 is a rapidly 
rotating, hot, massive, hot white dwarf, formed in the 
merger of two normal white dwarfs. Such a rapidly rotat- 
ing white dwarf can have sufficient rotational energy loss 
(~ 10 36 ergs _1 ) to power the observed X-ray flux. Using 
blackbody emission and a white dwarf radius of 0.002 Rq, 
our optical limits correspond to a maximum temperature 
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Fig. 1. R and I band images around IE 2259+58.6, the circle is the 95% confidence circle obtained in Sect. |[ This circle 
is also shown in the lower three enlarged images, where one has the star identifications and two have the conservative 
Rosat error circle at 95% confidence and previous error circles obtained by Fahhnan et al. (1982; short dashed) and 
Hanson et al. (198S; long dashed) superposed. 



kT = 0.4 keV. Since the Eddington c ondition gives a max- 
imum temperature of kT = 0.2 keV ( Paczynski 1990 ), our 



tion, due to the irradiation. The second model we use is 



data c; ji thus not exclude this model. 



Similarly, our limits do not allow us to constrain the 
magnetar model. (A very crude estimate gives R ~ 37 for 
a kT — 0.4 keV blackbody at d — 6 kpc and reddened by 
A R = 3.8.) 

5.2. An isolated NS accreting from a disk 

We will make an estimate for the optical emission of an 
accretion disk using two simple models. In the first, we use 
the T{r) oc r~ 3 / 7 relation expected for an accretion disk 



dominated by irradiation (labeled V90 hereafter; Vrtilek 



et al. 199C ). An important assumption in obtaining this 
relation is that the disk is isothermal in the vertical direc- 



a Shakura-Sunyaev disk (labeled SS; 3hakura & Sunyaev 



1973), in which the energy generated by the disk itself is 



taken into account but irradiation is neglected. For both 
models, we assume the disk emits like a blackbody. We 
believe the first model should be the more reasonable ap- 
proximation (at least for radii larger than a few 10 10 cm; 
for smaller radii, this model underestimates the emission; 
see Fig. ||), while the second model should give a conser- 
vative limit. 

Fig. H shows predicted R-band fluxes as a function of 
disk radius for a disk inclined by 60° relative to the line of 
sight, located at a distance of 6 kpc, and either irradiated 
by a central X-ray source with Lx = 4.3 x 10 35 ergs _1 
(V90), or not irradiated (SS). Also shown is the contri- 
bution to the total flux from each radius r (at constant 
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bg(disk radius (cm)) 

Fig. 2. Cumulative R band flux observed at earth as a 
function of radius for irradiated (solid line) and unirradi- 
ated (dashed line) disks at a distance of 6 kpc. The lower 
curves show the contribution of the flux for individual disk 
annuli (at constant A log r) . No correction was made for 
reddening. 



Alogr). There is no significant contribution to the flux 
for radii larger than -10 12 cm (V90) and -10 n cm (SS), 
because the temperature becomes too low. For a large disk 
around an isolated neutron star, we find R = 22.7 for the 
V90 disk (r > 10 12 cm; i = 60°), and R = 24.7 for the SS 
disk (r > 10 11 cm; i = 60°). Here, for both estimates we 
used an R-band extinction Ar = 3.8, which corresponds 



to Ay = 4.7 (see Sect. 2.1) 



From the above limits, we conclude that it is unlikely 
that IE 2259+58.6 is an isolated NS accreting from a disk. 
Only if we underestimated the reddening considerably or 
if the disk inclination were very unfavourable, a (large) 
accretion disk could have escaped detection. Increasing the 
distance to the source is not an option, since the optical 
emission of the disk is proportional to the X-ray emission, 
(for the SS disk model at lOkpc, we still expect R = 25.1). 

5.3. A very low mass binary 

The optical emission of low mass X-ray binaries is domi- 
nated by the emission from the accretion disk surround- 
ing the compact object. An empirical relationship, My = 
(1.57+0.24) - (2.27+0.32) x log[(L x /£cdd) 1/2 (i7l hr) 2 / 3 l- 



was found by van Paradijs & McClintock ( 1994| ). Assum- 
ing (V — R)o ~ 0, according to this relationship a puta- 
tive binary must have P ^ 10 2 s in order to remain below 
our detection limit R > 25.7. This would rule out all but 
a white dwarf counterpart. There is considerable scatter 
around the relation, however, and period estimates will be 
off by a factor of 5 for a change of 1 mag. Furthermore, it 



is far from clear one can extrapolate their result to such 
short periods. 

Given the above, it is perhaps more relevant to com- 
pare our limits with the emission from the compact low- 
mass X-ray binary 4U 1626-67 ( |Chakrabarty 1 998|). In- 



deed, this source has a num b er of the properties of AXP 
QMcrcghctti fc Stella 1995] cf. |van Paradijs ct al. 19951 ). 
Under the assumption that 4U 1626-67 and IE 2259+58.6 
have the same Lx/L opt ratio, the expected magnitude 
would be R = 24.8. Here, we have taken account of the dif- 
ference in X-ray flux of about a factor ten (compare, e.g., 
|Orlandini et al. 199£ and Parmar et al. 1998| ), as well as 
the difference in reddening towards the two sources. Even 
if the X-ray flux of IE 2259+58.6 were only one twentieth 
that of 4U 1626—67, we would still expect to have detected 
the optical counterpart at R = 25.5. Thus, we conclude 
that most likely even such a compact binary would have 
been detected. 

What would our limits imply if there were no accre- 
tion disk present and the optical emission was dominated 
by the companion star? This seems unlikely, but might 
perhaps be the case if accretion is from a (irradiation- 
induced) stellar wind. We will discuss two possible com- 
panions below: (i) a main-sequence star; (ii) a He-burning 
star. From Sect. 5.1) above, it will be clear that we can 
put no meaningful constraints for the case the only source 
of optical emission is a white-dwarf companion. 

Based on the colors derived from model spectra by 
Bessell et al. (199S) (overshoot models, their Table 1, for 
logy = 4.5), we find that a main-sequence companion to 
IE 2259+58.6 would have to be of type K or later to escape 
detection. 

To start helium burning, a star's mass must exceed 
^0.3Mq. From Savonije et al. ( 1986 ), we find that a ~ 
O.3M helium burning star has T cS ~ 2 0000 K and L ~ 
0.6L©. Using Table 3 of Bessell et al. (|1998| ), we infer 
a bolometric correction of —2.6 and V — R = —0.1 for 
a 20000 K star (note that this table does not extend to 
logg > 5, but since the variation with logy is small, our 
extrapolation should be safe). For these values, the star 
would have R = 25.7 and / = 24.7, which is just below 
our detection limit. More massive stars would be brighter, 
and, thus, if the companion is a star burning helium it 
must have a mass close to the limit of 0.3 Mq. 



6. Conclusions 

We have shown that the tight limits we derive on the op- 
tical emission of IE 2259+58.6 are a serious problem for 
scenarios involving accretion by way of an accretion disk. 
This is irrespective of whether the source is in a binary, 
except if the binary is extremely compact. Therefore, we 
conclude that, most likely, IE 2259+58.6 is an isolated ob- 
ject, and that its emission is not powered by accretion. We 
note that the same conclusion was indepently drawn by 
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Kaspi et al. (1999) on the basis of the extreme rotational 
stability of IE 2259+58.6. 
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